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Abstract 

We elaborate the idea that the matrix models equipped with the gauge symmetry provide a 
natural framework to describe identical particles. After demonstrating the general prescrip- 
tion, we study an exactly solvable harmonic oscillator type gauged matrix model. The model 
gives a generalization of the Calogero-Sutherland system where the strength of the inverse 
square potential is not fixed but dynamical bounded by below. 
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1 Introduction 



Recent advances in the noncommutative field theory have enabled us to realize that the quantum 
Hall system is closely related to the noncommutative Chern-Simons theories. After the pioneering 
work by Susskind [Q], utilizing the fact that the noncommutative field theories can be formulated 
by matrices of infinite size, Polychronakos proposed a finite matrix Chern-Simons model for the 
description of the finite number of electrons in the quantum Hall system M. Soon after, a matrix 
version of the Laughlin's wavefunction [[| was recovered by Hellerman and Raamsdonk [Q]. 

In this paper, we wish to clarify the underlying principles of the above matrix model 
approach to the quantum Hall system in a more general setup. We elaborate the idea that the 
matrix models equipped with gauge symmetry provide a natural framework to describe identical 
particles (see e.g. || ||). 

After demonstrating the general prescription, we explicitly study an exactly solvable har- 
monic oscillator type gauged matrix model. We show the model gives a generalization of the 
Calogero-Sutherland system where the strength of the inverse square potential is not fixed but 
dynamical bounded by below. 



2 Description of Identical Particles 



One of the intrinsic properties of the fundamental particles in nature is the very fact that they 
are identical. Namely it is in principle impossible to identify each individual particle at different 
time slices. In the ordinary quantum mechanics, the conventional way of incorporating this idea 
is to anti-symmetrize the wavefunctions over the particle indices by hand. 

Matrix gauge theories provide a more natural framework to describe identical particles. To 
see this we first consider the classical Lagrangian for N electrons 

N 

L = Y^ T(x a , x a ) - V(x 1 ,x 2 , ■ ■ ■ , x N ) , (1) 

a=l 

where x G R/°, the D-dimensional 'space', and the potential is totally symmetric over the particle 
indices. Formally introducing a diagonal N x N matrix, X = diag(xi, x 2 , • • • , xn), one can always 
rewrite the Lagrangian in terms of the matrix 



tr 



T(X,X) -V(tr/,(X)), (2) 



where we used the fact that the symmetric potential can be always written in terms of the 
traces of certain set of functions, tr//(X) = J2a=i fi( x a)- We put the proof at the end of this 
section. In particular, when x a carries no spatial index as in the D = 1 case, one can simply set 
flip) = x l , 1 < I < N. 

Allowing the off-diagonal elements with the restriction, X being hermitian to ensure the 
eigenvalues to be real, we encounter a new physical system. First, the action acquires a U(iV) 
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symmetry for which the matrices transform in the adjoint representation 

x -> uxW , u e u(JV) . (3) 

We note that the permutation group relevant to the relabeling of the electrons is a subgroup of 
\J(N). As electrons are identical, the labeling is not physical so that the permutation symmetry 
is an auxiliary one. This suggests to gauge away the U(iV) symmetry in the matrix formalism. 
Namely we introduce an auxiliary matrix, A and replace the ordinary time derivative by a 
covariant time derivative 

D t X = X + i[A ,X]. (4) 

Then what A brings new is the Gauss' constraint or the equation of the motion for A , which 
gives the quantum generator of the U(iV) symmetry from the Noether theorem. Since we embed 
the discrete permutation group into the continuous group, the reparameterization of the electrons 
can be now realized by quantum operators. At the quantum level the constraint is to be imposed 
on the wavefunctions and this will mode out the auxiliary symmetry. We emphasize the point 
that the U(A r ) symmetry (||) can be now time dependent, and physically this amounts that we are 
requiring the physics to be invariant under not only time independent but also time dependent 
reparameterization of the particles. Once again, it is in principle impossible or meaningless to 
identify each individual particle at different time slices. 

One interesting "generalization" is to add a term into the action which is linear in A [||, |?J 

K barc tr^4 . (5) 

While writing the Gauss' constraint at the quantum level there always occurs an ordering ambigu- 
ity, since the constraint contains the product of X and its conjugate momentum. The ambiguity 
amounts to adding an identity matrix to the Gauss' constraint up to a factor. Thus, considering 
the above term (||) is not a mere "generalization" but rather a natural requirement. After all, the 
final form of the Gauss' constraint which is to be imposed on the quantum wavefunctions should 
be written in the normal ordered way with the physical coefficient, k 

U(iV) generator = kI nxn . (6) 

" V ' 

normal ordered 

After all, the actual integer value of the physical coefficient, k is given by the superselction rule. 

The consistency at the quantum level requires the physical coefficient, k to be an integer 
H 0- At the quantum level, the Gauss' constraint generates unitary transformations, U = 
e lA , At = A on all the arguments in the wavefunction [Q] 

|*') = e iKtrA |$). (7) 

Taking the particular choice, A = diag(2-7r, 0, 0, • • • , 0) gives the identity matrix, U = 1 NX n, 
and apparently the Gauss' constraint on the wavefunctions works successfully only for integer, 
k. Essentially this quantization is identical to that of the coefficient in the noncommutative 
Chern-Simons theories ||. 

Nevertheless, non-zero k is yet problematic for the finite matrix models. As the matrices 
are in the adjoint representation, the central U(l) transformation would leave the wavefunction 
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invariant, and this is clearly inconsistent with Eq.(0) for non-zero k. Curing the problem requires 
the presence of new variables in other representations. A natural candidate is a complex bosonic 
vector, (j> in the fundamental representation so that D t (f) = eft + iA (ft. As we will see later, this 
new variable governs the strength of the intrinsically existing repulsive potential in the matrix 
gauge theories. 

Introducing the complex vector, the gauged matrix model can be self-consistent and is now 
of the general form 

L = tr[T(X, D t X) + K baro A ] - V (tr/ a (X)) + vector parts . (8) 

As worked out in [|], [?J the general quantum wavefunction satisfying the Gauss' constraint consists 
of U(iV) invariant part times |k| products of SU(A r ) invariant parts (see Eq.(|l9|)). Due to the 
latter, the wavefunction is an eigenstate of the particle exchange operator of the eigenvalue, (— 1) K 
capturing the identical nature of the particles. 

The prescription to give the physical meaning to any expectation value in the matrix model 
is the map 

(tt|tr/(X)|tf) «— Jdx D p(x)f(x) , (9) 

where p(x) is the corresponding physical density function, while f(x) is an arbitrary function. 
As we turn on the off-diagonal components in X and introduce the vectors, the present non- 
matrix system is not simply equal to the initial one (|]). Physically, the off-diagonal element, 
X a b, a 7^ b corresponds to the interaction between the two particles, a, b, and integrating out the 
off-diagonal components generates an inverse square type repulsive potential. Thus the matrix 
gauge theories intrinsically contain an repulsive potential among the particles. However, this can 
be in principle eliminated in the matrix model by adding the counter term written in the matrix 
form, if necessary. 

It is worth to note that due to the Gauss' constraint on the wavefunction (|^) any matrix 
valued expectation must be U(A^) invariant, and hence 

(*| f(X) |*) = | tr/(X) |*) I NXN . (10) 

For the systems of kinetic terms quadratic in time derivatives in the D dimensional space the 
phase space for the matrices has the dimension, 2DN 2 so that subtracting the Gauss' constraint 
and the gauge symmetry, the physical degrees of freedom for the matrices is 2{D — 1)N 2 + 1. 
On the other hand the vector has the linear in N degree of freedom. Hence when D = 1 the 
total degrees of freedom is linear in N, otherwise it is quadratic in N. Nevertheless if we include 
the BFSS type potential, tr[X / ,X J ] 2 [§, at the low energy limit the matrices tend to commute 
each other making the "off-diagonal components" negligible. This will result in linear degrees of 
freedom at low energy. In fact, in a special low energy limit the particles form a collinear motion 
j?], [lOfl . On the other hand, for the kinetic term linear in time derivatives like the D = 2 matrix 
Chern-Simons model, the dimension of the phase space is halved and the total degrees of freedom 
in linear in N. 
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The proof of the fact 

Consider a symmetric function depending on N particles' coordinates, x a ,a = 1,2,---,JV of the 
form 

F(xi,---,x N ) =^2 fi{x pi )f 2 (x P2 ) ■ ■ ■ f N (x PN ) , (11) 
p 

where /i, ■ ■ ■ , f N are functions of one particle coordinate and the sum is over the N\ permuta- 
tions so that the function is apparently symmetric over the particle indices. General symmetric 
functions can be written in terms of this kind of symmetric functions. For example, the Coulomb 
interaction can be written as a fraction of such two functions. In the below we show that F can 
be written in terms of tr fi(X) with X = diag(xi, x 2 , ■ ■ ■ , x N ). We prove this by the mathematical 
induction on the number of non-constant functions in / 2 , • • • , f N } which we denote by If 
#f = 1, it is easy to see 

F(x u - ■ ■ ,x N ) = (N - iy.tr f a (X)H f p , (12) 

where f a is the only one non-constant function. Hence the statement holds for #j = 1. Now 
we assume that F can be written in terms of tr/;(X) for #j < n cases, and consider F in the 
#/ = n case. We let with out loss of generality fi,-- - ,f n be the n non-constant functions and 
set 

F< = F-{N-n)\ (f[tvf a (X)) [Hfp). (13) 

\a=l ) \/3>n J 

It is crucial to note that F 1 belongs to the classes, #/<nso that F' and hence F can be written 
in terms of tr fi(X). This completes our proof. 



3 Generalization of Calogero- Sutherland System 



The D = 1 gauged matrix model we consider here is the harmonic oscillator type. With the 
column vector, <fi, and the row vector, <ft = (jr , the Lagrangian is 



tr 



\m{D t Xf - \gX 2 + mD t <j>D t <i> - g' <M> + ^ barc A . (14) 



If g = g', writing a (N + 1) x (N + 1) matrix {ft), the model can be regarded as the truncation 
of the bigger matrix model with the broken gauge symmetry, ( „ ° ), \J(N + 1) — ► U(iV). 

With the conjugate momenta, P = mD t X, = mD t cj) we define the following dimension- 
less quantities 

C = -j= [{mg)i.X + i{mg)'\p^j , rj± = ^= ((mg')*(f) ± i{mg')-^p^j . (15) 

The standard quantization shows the non-vanishing commutators are 

[C\,C c d ] =5 a d 5 b c , [r) a + ,fj +b ]=5 a b , [ V a ^,fj^ b } = -5 a b . (16) 

Thus, C, fj + ,rf- and C, rj + ,fj- are respectively creation and annihilation operators (our somewhat 
unconventional notation for r/_ is to keep the consistent U(iV) index notation.). The appearance 
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of the two sets of vector-valued harmonic oscillators is due to the fact that the kinetic term for 
the vector is written with the quadratic time derivatives. By introducing one more vector field, 
one could equivalently rewrite the vector parts in the action to contain two kinetic terms linear in 
time derivatives having the opposite signs (see [11] for further generalization of the kinetic terms 
linear in time derivatives). 



In terms of the operators, the Gauss' constraint reads in the normal ordered form, 

c\c\ - c\c c h + n+b v a + - v- n-b = ^s a b , 



(17) 



where the left hand side generates the U(iV) transformations, while k on the right hand side is 
now of the physical value. The Hamiltonian is 



H 



tr 



{CC + CC) + Ji{r ]+ f) + +ri-n- 



(18) 



£tr(<7C) + (77+77+ + tr(77„77_)) + 



3. N 2 + 

m 



3-N. 

m 



Here the second expression is written in the normal ordered fashion so that tr(CC), 77+77+, 

tr(77_77_) are the number operators and \\J^N 2 + ylj^N is the zero-point fluctuation of the 
energy. 



The exact wavefunctions satisfying the U(JV) covariance condition ([7|) due to the Gauss' 
constraint are for k > and k < respectively M 

f Gx (e a ^" a »fj +ai {f)+C) a2 ■ ■ ■ {f) + C N - l ) aN ) R |0) 

I*) = (19) 

where G = G{trC l ,f] + C m 77_) is an arbitrary function of the U(iV) invariant building blocks which 
are made of the creation operators only, tvC 1 and fj + C m r]- with 1 < I < N, < m. The trivial 
G corresponds to the ground state whose energy is 

£0 = \ -(M + 1)N + I x [^((\k\ + 1)N 2 - \k\N) . (20) 

In particular when n = 0, the vacuum, |0) is the ground state, and in this case we can 
calculate the exact density function, po(x) in Eq.(^). First, using the large N behaviour, 

(0|tr(C + C) 2m \0) ~ (2m ] ! N m+1 , (21) 
m\(m + lj! 

one can check that the Fourier mode for po(x) is given by the Bessel function, J\. Realizing the 
Fourier transformation of a half circle is J\ we get for large N 



7T 



Po ( x ) = lVw(mg) 1 iRe(Jl-^xA . (22) 
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In fact, this density function is identical to that in the Calogero-Sutherland system [12]. An 
intuitive way to see this result is to note the close relation to the matrix quantum Hall system 
with the confining harmonic potential [0], where the density function is constant on a disc. 
Essentially our density function is the one-dimensional projection of it. For general k at the 
center, x = we expect 

hm hML= f* (23 ) 
N^oo 2N tt 2 (|k| + 1)' V ; 

and this is an analogue of the fractional filling factor, v = 1/(|k| + 1) in the quantum Hall system. 
The "+1" in the denominator is again due to the zero-point fluctuation or the Vandermonde 
determinant. 

Henceforth we discuss the classical dynamics of the system focusing on the k > case. The 
equations of motion are 

= mD t D t X + gX , (24) 
= mD t D t cj) + g'4> ( O = D tV± ± iJ^r]±) ■ (25) 



As in Q we choose the gauge such that X is diagonal and 77* is real and non-negative. For 
the simplicity of notation we define a matrix, K a b = r/" r/+b — rft_f]_b ■ Now the Gauss' constraint 
determines r]^_ and the off-diagonal components of A completely 

K a b /m 



tI% = Jk+\t£\2, A a 0b = - ^— for a + b. (26) 

Substituting these into the Hamiltonian we obtain a generalized Calogero-Sutherland model 
(cf.0) 

= E + h-l + ^n-av a - + E » (27) 

o=l a>b ^ a Jbb) 

where K a b is now a function of r/_ only 



K a b = ^J(K+\V-\ 2 ){K + \v b -\ 2 )- V-V-b ■ (28) 

Apparently the complex vector, r/^. gives the novelty. Namely the strength of the repulsive 
potential, \K a b\ between two particles is not fixed but varies. The Schwarz inequality shows 
\K a b\ > K an d the saturation occurs when rft_ is independent of the particle index, a. In partic- 
ular, for the finite energy configurations with the large g' limit, r/_ vanishes classically and the 
saturation occurs. 

The corresponding Lagrangian is with (p£|) 

X — \K a b \ 2 /m 

a=l a>b 



\Xa Xb ) 
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The minus sign for the kinetic term of r/_ is essentially from Eq.(|iq). Consistency requires that 
the dynamics of this generalized Calogero-Sutherland model must agree with the full equations of 



motion of the matrix model ( |24| , 25). Some straightforward manipulation with Eq.(^fj) can show 



that only when the diagonal component of A is given by 

A U = ~ £ U A > + A b 0a )(n b + /r,l) , (30) 

we get Dtrj + + i\J~^rj+ = 0, and the matrix equations of motion reduce to those of the present 
generalized Calogero-Sutherland model 

_ 2\K\\ 2 /m 
= mx a + gx a - -3- , 

b+a [ - Xa Xb > 



(31) 

— \K a A 2 I™ 

o = ^ + t£ + £ 



wr a \K a b\ 2 /™ 



b^a ( Xa Xf> ) 

The former comes from the diagonal components in (^4|) while the off-diagonal components vanish 
due to (|D, (mD t D t X + gX) a b = -i(D t K)\/(x a -x b )=0,a^ b. 



^_=ce 2 ^- {t+U \ (32) 



Especially when the effective charge saturates, \K a b \ = k, we get the solution for the two- 
particle system in a closed form 



and xi, X2 are 



It would be interesting to see any quantum correction to the inverse square potential in 
Eq.(29). When the vector freezes it should reduce to k 2 — » k(k + 1) In any case, since k is 



l CM sin(, /£ t) ± iJsmpJ*- 1 + 0) + ^ + k?/(4gmP) . (33) 



an integer, the charge of the potential is quantized. 



4 Conclusion 



We have argued that the gauged matrix models provide a natural framework to describe identical 
particles since the particle indices which are not physical are gauged. The models intrinsically 
contain the interaction among the identical particles. Especially for D = 1 or the [X^X* 7 ] — > 
limit, it corresponds to the inverse square type repulsive potential. As an example we have 
considered an exactly solvable harmonic oscillator type gauged matrix model. We have obtained 
all the exact quantum wavefunctions and demonstrated that the model is a generalization of the 
Calogero-Sutherland system where the strength or the charge of the potential is not fixed but 
dynamical bounded by below. Furthermore, when the saturation occurs the charge is quantized. 

The model we considered is non-relativistic as the kinetic term is quadratic in time deriva- 
tives, and this is essentially the reason why we get the inverse square potential. It is very much 
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desirable to find a matrix gauge model where the intrinsic interaction is given by the inverse 
potential or the Coulomb interaction. In this sense, it is worth to study the relativistic matrix 
model and analyze its intrinsic interaction such as 

S = J dTtTyf(D T T) 2 - (D T X) 2 . (34) 

Quantization of the system is analogue to that of the string theory. It would be interesting to see 
if the critical dimensions free of the Lorentz symmetry anomaly exists. A single particle system 
i.e. lxl matrix model is anomaly free in any dimension. On the other hand, the generic matrix 
model including the off diagonal elements may be anomalous just like the string theories where 
the anomaly comes from the oscillators not from the zero mode. Eq.(|34|) is a straightforward 
generalization of the single relativistic particle action, and an alternative approach to write the 
relativistic matrix model would be the dimensional reduction of the non-Abelian Born-Infeld 



action fl5|| . 



An electron in one place seems to be "distinguishable" from one at far distance. To explain 
this large scale phenomenon, the spontaneous symmetry breaking of the U(iV). is to be considered. 
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